Interplay between non equilibrium and equilibrium spin torque 
using synthetic ferrimagnets 
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We discuss the current induced magnetization dynamics of spin valves Fo|N|SyF where the free 
layer is a synthetic ferrimagnet SyF made of two ferromagnetic layers Fi and F2 coupled by RKKY 
exchange coupling. In the interesting situation where the magnetic moment of the outer layer 
F2 dominates the magnetization of the ferrimagnet, we find that the sign of the effective spin 
torque exerted on the free middle layer Fi is controlled by the strength of the RKKY coupling: 
for weak coupling one recovers the usual situation where spin torque tends to, say, anti-align the 
magnetization of Fi with respect to the pinned layer Fo. However for large coupling the situation 
is reversed and the spin torque tends to align Fi with respect to Fo. Careful numerical simulations 
in the intermediate coupling regime reveal that the competition between these two incompatible 
limits leads generically to spin torque oscillator (STO) behavior. The STO is found in the absence 
of magnetic field, with very significant amplitude of oscillations and frequencies up to 50 GHz or 
higher. 

PACS numbers: 72.25.Ba, 75.47.-m, 75.70.Cn, 85.75.-d 



Since the first successful experimental manipulation of 
magnetic configurations by spin-polarized currents [l|, Q , 
the interest in spintronics devices entirely controlle d by 
electrical currents has rapidly increased [§-flol \vX fl3j . 
The key concept behind theses experiments is the no- 
tion of spin transfer torque (STT) predicted by Slon- 
czewski |14| and Berger [15| : a current gets spin polarized 
by a first (pinned) magnetic layer and then transfers some 
of its polarization to a second free layer, hence exerting 
a torque on its magnetization. As a result the magneti- 
zation can switch to a different static configuration 0, Q 
or even to a dynamical stationary regime where the mag- 
netization of the free layers shows a sustained preces- 
sion [H-Q- This dynamical regime, known as the spin 
torque oscillation (STO), allows to convert a dc voltage 
into a microwaves signal at several GHz and has a real 
potential for applications (nanoscale microwave source 
or detector, high frequencies possibly above CMOS tech- 
nology, narrow band). There are still technological is- 
sues with STO based devices, however, such as very lim- 
ited output power and the need of rather strong external 
fields. Indeed, the main mechanism for producing STO 
is based on a detailed tuning of the angular dependancc 
of the spin torque [B| and requires the use of an exter- 
nal magnetic field. Recently, different mechanisms have 
emerged to enable STO behavior in the absence of exter- 
nal field: the first approach [§] relies on the precession 
of a non uniform magnetic structure, a magnetic vortex, 
that can be excited at sub GHZ frequencies; a second 
approach [H [n| uses a strongly asymmetric spin valve to 
increase the anharmonicity of the angular dependance of 
the spin torque, leading to the so-called "wavy" behavior 



of the ST T ll6H19j . Other approaches use an orthogonal 
polarizer jlCf or a perpendicular free layer [ll. 12|. 

In this letter we propose an entire new mechanism 
which naturally induces strong STO behavior, even in 
the absence of magnetic field, using a synthetic ferrimag- 
net (SyF) as the free layer of the spin valve. The SyF is 
made of two ferromagnetic layers coupled thr oug h anti- 
ferromagnetic RKKY exchange interaction Here, 
we predict that by simply tuning the effective strength 
of the RKKY exchange coupling (through the different 
thicknesses of the layers for instance), the sign of the 
effective STT can be changed. At intermediate effective 
coupling, the STT generically induces strong STO be- 
havior irrespectively of the presence of applied magnetic 
field. 

In the following we first provide a simple physical pic- 
ture that captures the essential features of our setup. 
Our main prediction uses very general conservation ar- 
guments and should therefore be quite robust and inde- 
pendent of the details of the model. In a second step, 
we provide more quantitative arguments to describe the 
strong coupling regime. Last, we present a detailed nu- 
merical study of the phase diagram of our model. We 
find strong STO behavior at rather large frequencies in a 
wide portion of the phase diagram hinting that SyF based 
STOs could be good candidates for application purposes. 

Physical mechanism. We consider a nanopillar spin 
valve Fo I No I SyF made of a pinned ferromagnetic layer Fo 
and a free synthetic ferrimagnetic layer SyF = Fi|Ni|F2 
(where No and Ni are normal spacers), see Fig. QJi for a 
cartoon of the full stack. Upon injecting a current den- 
sity j through the nanopillar, a spin torque r is exerted 
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Figure 1: Panel (a) is a cartoon of the magnetic multilayer 
stack, in which Mi is the the magnetization of the layer Fi. 
The angle between Mi and Mo (Ma) is respectively 6 (a). 
Panels (b,c,d) show the stationary magnetization mf along 
Mo as a function of the current density j in A.cm~ 2 , for a 
stack with di = 2. nm, c?2 = 4 d\ and three different coupling 
strength J; (b) weak coupling limit (J = 0), (c) Intermediate 
coupling limit (J = — 10~ 3 J.m~ 2 ) (d) Strong coupling limit 
(J = —6. 10 -3 J.m -2 ) Upper (Lower) inset : Sketch of the 
torque Ti action on Mi (Ms y f) for the weak (strong) coupling 
regime. 



on the magnetic moments of the layers |23( . This torque 
has been extensively discussed in the literature and can 
be strong enough to destabilize the initial magnetic con- 
figuration. On the other hand, the magnetization Mi 
and M2 of Fi and F2 are coupled by an oscillatory ex- 
change (RKKY) interaction J which is nothing but the 
spin torque present in the pillar at equilibrium [/Ml . |25| . 
The coupling J is reminiscent of Fricdel oscillations and 
typically behaves as J ~ (cos2kpdN)/dN where &f is 
the Fermi momentum and cIn the width of the spacer 
Ni. J can be tuned from negative (antifcrromagnctic 
coupling) to positive (ferromagnetic coupling), hence the 
name " oscillatory" . Here we focus on negative values of 
J which stabilize an anti-aligncd configuration of Mi and 
M 2 . The goal of this letter is to discuss new phenomena 
that arise from the competition between the equilibrium 
and the non equilibrium torques. 

In order to gain a qualitative understanding of the un- 
derlying physics, let us discuss two extreme limits. Wc 
start from a configuration where Mi is aligned with Mo 



while M2 is anti-aligned, and inject electrons from the 
right. Upon crossing the pinned layer Fo, the electrons 
get a polarization anti-parallel to Mo. Denoting J; the 
spin current (per unit surface) flowing in the normal layer 
N i; we have Jo = hpjm /(2e) with m = M /|M | and 
p the polarization of the current (0 < p < 1). Let us fo- 
cus on the weak coupling limit where J — > 0. When there 
is a finite angle 9 between Mo and Mi, the spin current 
on the left Jo and right Ji of Fi become different and 
a torque T\ = Jq — J± is exerted on the magnetization 
of Fi: conservation of angular momentum implies that 
whatever spin lost by the conducting electrons is gained 
by magnetic degrees of freedom (l4| . This conservation 
of angular momentum reads, 
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Tl 
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where d\ is the width of Fi, 71 its gyromagnctic ratio 
and the . . . indicate other contribution to the dynam- 
ics (magnetic anisotropy, damping) to be discussed later 
(upper inset of Fig. [1]). Ignoring (momentarily [26j]) the 
role of F 2 , we arrive at 88 /dt = jihpj sin#/(2eeJi|Mi|) 
and recover the usual phenomenology of spin torque: the 
torque tends to stabilize the configuration where Mi is 
anti-parallel to Mo (or destabilize it if one reverses the 
sign of the current). 

Let us now discuss how this picture is modified when 
one considers the opposite strong coupling J — > — 00 
limit. Upon switching on J, the exchange energy (per 
unit surface of the pillar) Ej = — J (mi • m 2 ) induces 
a field like term in the RHS of Eq.(Q]) of the form 
+ Jm 2 x mi. One needs also to consider the correspond- 
ing equation for the dynamics of M 2 (see below Eq.([3]) 
for the full model) and one obtains a potentially rich cou- 
pled dynamics for Mi and M 2 . The situation simplifies 
in the limit where the exchange energy J — > — 00 domi- 
nates: the relative configuration of Mi and M 2 becomes 
frozen in the anti-aligned position. Conservation of an- 
gular momentum now leads to a spin torque on the total 
magnetic moment: 



d_ 

dt 



diMi 

7i 



d 2 M 2 



")2 



(2) 



where we have used the fact that r 2 = Ji — J 2 vanishes 
when Mi and M 2 arc collincar. Note that the exchange 
field, which conserves the total angular momentum, re- 
distributes the torque Ti on Mi and M 2 so that they 
remain anti-aligned. 

Eq.([2]) allows two very different regimes: if 
di|Mi|/7i > d 2 |M 2 |/7 2 then the effective magnetiza- 
tion direction Ms y F of the SyF layer points along Mi 
and the torque is very similar to the weak coupling 
regime. This case has attracted some attention recently 
both experimentally by Smith et al. (2?J and theoreti- 
cally by Balaz et al. [28[. Here, we focus on the other 
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limit di | Mi I/71 < d2\M-2\/l2 where Ms y F points in the 
opposite direction as Mi (lower inset of Fig. [I}. In 
this limit, the torque tends to stabilize the configura- 
tion where Mg y p is anti-aligned with M hence where 
Mo and Mi are aligned: this is the opposite situation 
from the weak coupling J — > limit. 

To summarize, when going from J — > to J — > —00, 
without changing the sign of the current, the torque tends 
to favor two different stationary situations. In the rest 
of this letter, we study this crossover in more details. In 
particular, we find that the intermediate coupling regime 
reveals interesting, STO like, dynamical behavior. 

Model. We now turn to the modclisation of our sys- 
tem. The full stack we arc considering has the form (see 
Fig. UK): AF|F |N |Fi|Ni|F 2 |N2, where AF is an antifer- 
romagnetic layer (typically IrMn) that pins the magne- 
tization of Fo . We set our reference spin axis e z parallel 
to Mo. Magnetization dynamics is described by two cou- 
pled Landau Lifshiz Gilbert (LLG) equations that read 
(in SI units), 

dmi dmi 71 . . 

: 71B1 x mi + o-imi x — — + t 1 (3a) 



dt 
dm 2 



7 2 B 2 x m 2 + a 2 m 2 x 



dt |Mi|di 
dm 2 72 
dt \M 2 \d 2 



r 2 ,(3b) 



where nij = Mj/|Mj| is the unit vector describing the 
magnetization orientation of the layer Fi and a>i is the 
damping factor of the corresponding magnetic material. 
The effective field Bi seen by the magnetization in the 
layer (i = (1, 2)) is 
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where i = 3 — i corresponds to the index of the other 
layer. The effective field includes an uniaxial anisotropy 
K % u field and the RKKY exchange field. In the physical 
regimes studied below, RRKY coupling is the dominat- 
ing energy so that our results are rather insensitive to 
the presence of less important terms (for instance dipo- 
lar coupling between layers that we do not take into ac- 
count). The spin torque Ti is calculated in the framework 
of CRMT (Continuous Random Matrix Theory) [HEH. 
CRMT is a semiclassical approach that generalizes Valet- 
Fert theory to non collinear situations. For discrete sys- 
tems, It is formally equivalent to the generalized circuit 
theory [3(| and is well adapted to the treatment of metal- 
lic magnetic multilayers. 

Strong coupling limit. When the RKKY coupling is 
large, one can derive an effective LLG equation for the 
dynamics of the SyF. It reads, 



<9mi 
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~~dT 



7cff 



-n, 



(5) 



where the various effective parameters get renormalized 
as follows: %s = 7i/(l " 8), "off = ("i + 5a 2 )/(l - S), 



Kf = Kl + K 2 u d 2 /d u B cff = (JCf/lMiDImreJe, and 
the rcnormalization parameter 5 has been defined as, 



71IM2M2 
72 Mi di 



(6) 



(Eq.© is obtained by calculating d/dt\ta.\ + <5m2] which 
cancels the RKKY contribution, and then setting 1112 = 
—mi and r 2 = 0). From this equation we can clearly 
see the inversion of the sign of the effective torque 7 e ffTi 
discussed previously when 5 > 1. 
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Figure 2: Upper: phase diagram of the system as a func- 
tion of the current density j and the thickness di for a fixed 
5 = 4 and J = —5. 1CP' ! J.m~ 2 . The various symbols indi- 
cate the observed stationary states in the corresponding re- 
gion (defined by different colors): t (l) stands for mf = 1 
(ml = —1), while STO corresponds to the presence a pre- 
cessional state. In the 2 STO region, two different STO 
states can be observed depending on the hysteresis history. 
Lower: Resistance R as function of time t for a pillar surface 
of 1000 nm 2 and d\ = 3nm. Upper (lower) plots correspond 
to initial conditions in the up (down) configuration. Left pan- 
els: j = 2.8 10 7 A.cm~ 2 . Right panels: j = 5.8 10 7 A.cm" 2 . 
The amplitude of the oscillation with the up initial condition 
is small and invisible on the scale of the graphics 



Numerical simulations. Let us now go back to the full 
model Eq.([3]) and study it numerically. We focus on a 
stack defined as IrMnslPyslCuiolCodi |Ruo.3|Cod 2 |Cuio 
where the index indicate the widths of the layers in 
nm. We have used K%° = 8. 10 4 J.m" 3 , M Co = 
1.42 10 6 A.m -1 and the transport parameters (spin re- 
solved bulk resistivities, interface resistivities, spin flip 
lengths) have been taken from the database established 
in MSU and CNRS/Thales laboratory [3l|. The coupling 
constant J can be tuned by changing the Ru thickness 
with typical values J Ru 



-5. 1CT 3 J.m" 2 [2l|. The 
torque and resistances have been calculated using CRMT 
and the corresponding parametrization (as a function of 
the angles 6 and a, see FigfTK) incorporated into our nu- 
merical LLG integrator. 

A first set of curve is presented in Fig. Q] b,c and d 
where the stationary magnetization m\ along the z axis 
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is plotted against the current density j for three values 
of the coupling constant: small (b), intermediate (c) and 
large coupling (d). For small coupling (b) , we recover 
the usual behavior of current induced magnetic reversal 
in spin valves: at zero current the system has two sta- 
ble configurations where mi is aligned (mf = 1) or anti 
aligned (mf = —1) with m . Upon injecting a strong 
positive current, one "pushes" toward the anti-aligned 
configuration which becomes stable above a critical value 
of the current density j (0.5 10 7 A.cm -2 in this example). 
At large coupling (d) , this hysteretic curve is reversed, in 
agreement with the arguments developed above. A very 
interesting regime is found at intermediate couplings (c) 
where one observes a rather large window where the sta- 
tionary value of mf corresponds to a finite angle between 
mo and mi , i.e. to a dynamical state where a sustained 
precession of mi is present (with frequencies around 15 
GHz in this example). 

The coupling J can be varied by tuning the width g?n 
of the normal (Ru) spacer. However, this is not very 
easy to control in practice, as the oscillatory character 
of the RKKY interaction make this variation non mono- 
tonic. Alternatively, one can explore the phase diagram 
of our stack by varying the thickness d\ for a fixed ra- 
tio 5 = dz/di, hence varying the relative importance of 
bulk versus surface terms in the LLG equation. Fig. [2j 
shows the resulting phase diagram in the (j, di) plane for 
S = 4. The different symbols indicate the possible sta- 
tionary states in the various regions of the phase diagram 
(defined by the various colors). The presence of two dif- 
ferent symbols (most regions) indicate that depending on 
the hysteresis history, one or the other state is observed. 
In particular, in the 2 STO region, two different STO 
states can be observed depending, for instance, on the ini- 
tial condition mf (t = 0) = ±1. These two different STO 
states merge upon entering the 1 STO region. We have 
investigated other stacks with different material parame- 
ters (not shown) and found very similar phase diagrams, 
including the presence of high frequency STO phases, as 
long as one remains in the S > 1 regime. We observe two 
kind of STO behaviors, as evident from the Resistance 
versus time R(t) traces plotted in Fig. [2j When the am- 
plitude of the oscillating signal is small (lower left of Fig. 
[5]), we observe some beating behaviors with one well de- 
fined frequency and a second much smaller frequency less 
well defined (corresponding to the precession of the two 
layers respectively) . When these two frequencies become 
closer, the time dependent signal gets bigger and R(t) 
becomes much more sinusoidal indicating phase locking 
between the precessional dynamics of the two magneti- 
zations mi and m.2- 

In the last Figj3l we focus on the region of the phase 
diagram (lower right corner of Fig|2]) where the ampli- 
tude of the STO signal is the highest. The upper curves 
show the amplitude of the time dependent signal while 
the lower curves show the corresponding frequency. We 



find that this large angle STO phase, which is stable in 
the absence of magnetic field, leads to very high frequen- 
cies up to several tens of GHz while sustaining high pre- 
cession angles (of the order of 7r/6 ). The observed high 
frequencies is tightly linked with the high value of the 
RRKY coupling which can be up to two orders of magni- 
tude larger than typical anisotropy or Zeeman energies. 
Even if the relative angle between the two magnetizations 
of the synthetic ferrimagnet remain close to n, this gives 
rise to frequencies of several tens of GHz. 

To conclude, we have shown that the interplay between 
equilibrium and non equilibrium torque in synthetic ferri- 
magnets based spin valves can lead to interesting physics 
and potentially a new route to high frequency STOs. 



< 



4.5 
3.4 
2.2 
1.1 










5:C 

: "if 
\ 


c 


( 










\ j [10* A/cm 2 


6 


0.0 





















| 5.6 














5 

j [10 T .4/cm- 


6.0 













1 2 3 4 5 6 7 



d 3 \nm] 

Figure 3: Upper plot: amplitude (in %o) of the time depen- 
dent signal of the resistance as a function of di for fixed 
j = 6. 10 7 A .cm . Lower plot: idem for the main oscillating 
frequency of the resistance. Insets: corresponding colorplots 
in the (j, di) plane. (5 = 4 and J — —5. 1CP 3 J.m -2 . 
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